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ABSTRACT 


A thermoelastic constitutive model for shape memory alloy (SMA) actuators and SMA hybrid 
composite (SMAHC) structures was recently implemented in the commercial finite element codes 
MSC.Nastran and ABAQUS. The model may be easily implemented in any code that has the 
capability for analysis of laminated composite structures with temperature dependent material 
properties. The model is also relatively easy to use and requires input of only fundamental 
engineering properties. A brief description of the model is presented, followed by discussion of 
implementation and usage in the commercial codes. Results are presented from static and dynamic 
analysis of SMAHC beams of two types; a beam clamped at each end and a cantilever beam. 
Nonlinear static (post-buckling) and random response analyses are demonstrated for the first 
specimen. Static deflection (shape) control is demonstrated for the cantilever beam. Approaches for 
modeling SMAHC material systems with embedded SMA in ribbon and small round wire product 
forms are demonstrated and compared. The results from the commercial codes are compared to 
those from a research code as validation of the commercial implementations; excellent correlation is 
achieved in all cases. 
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1. INTRODUCTION 

There has been growing interest in adaptive materials and structures in recent years. Structures that 
adapt to their environment or have the ability to change configurations have been of particular 
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interest. Shape memory alloys (SMA) have enormous potential for such applications. Their energy 
density is generally considered as the highest of known actuators, smart materials-based or 
conventional, so efficient actuation concepts may be possible. The interaction of temperature and 
stress on the transformation behavior of SMA materials can be used to exploit phenomena such as 
thermally induced shape memory, pseudoelasticity, as well as combinations of these effects. 
Significant progress has been made in the areas of SMA materials characterization, constitutive 
modeling, and integration into applications. However, structural concept development and 
component design have been stymied by the lack of an analytical tool in a general purpose structural 
analysis environment. Many structural concepts involving SMA materials have been developed by 
trial and error and/or using excessively simplified modeling approaches. 

There are various phenomenological and micromechanics-based constitutive models for SMA 
materials that are currently in use by various research groups . 1 ' 4 Although these models are very 
comprehensive in capturing the complex behavior of SMA materials, they are difficult to implement 
in commercial structural analysis codes and are difficult to use in practice. Thus, none of these 
models are in wide use, particularly by general structural analysts. An engineering model for SMA 
materials was recently developed based upon nonlinear thennoelasticity and definition of an 
effective coefficient of thennal expansion (ECTE ). 5 The ECTE model has the advantages of being 
relatively simple to use, readily implemented in commercial codes, and requiring input of only 
fundamental engineering properties. The ECTE model does, however, have a more restrictive range 
of application than the phenomenological/micromechanics-based models. Thus, it is envisioned that 
the ECTE model may be the model of choice for detailed design when appropriate and may be used 
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in conjunction with more a comprehensive model for analyses of increasing fidelity when the case 
demands. 

It is desirable in many applications to integrate SMA actuators into a structure, e.g., to achieve 
maximum space or weight efficiency, distributed or spatially varying actuation, autonomous 
response to environmental effects, minimum system complexity, etc. Previous studies have 
proposed special finite elements for modeling SMA materials according to various constitutive 
models, particularly the phenomenological/micromechanics-based models. These elements then 
provide the functionality of the actuator, separate from elements modeling the bulk of the structure. 
This approach is most effective when there are few actuators that are geometrically distinct from the 
distributed structure. This approach is inefficient, however, and intractable for practical structural 
analysis when there are many actuators, particularly when they are distributed within a host 
composite structure. This is especially true for cases involving SMA actuators of small size, like 
wire-type configurations. It is for these reasons that means for modeling SMAHC materials and 
formulating SMAHC elements is deemed necessary. The latter SMAHC modeling approach affords 
the ability to model an enormous number of material configurations efficiently, as long as the 
effective composite properties can be constructed from the constituents. The impetus for 
development of the ECTE model was the need to model SMAHC materials, so it is particularly 
suited to such modeling requirements. 

The concept of embedding SMA actuators in a composite structure was introduced by Rogers and 
Robershaw . 6 Several subsequent research efforts have investigated SMAHC laminates for a variety 
of applications . 7 " 1 1 The ECTE model has been used previously in a research code to analyze the 
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static and dynamic response of beam and plate structures. 5 ’ 12 ’^ It is the objective of this study to 
demonstrate SMAHC structural analysis using the ECTE model within the framework of 
commercial finite element codes. The thermally induced shape memory effect will be the focus of 
this work in application to structural thermoelastic, vibration, and shape control. A brief description 
of the modeling approach will be followed by discussion of implementation and use of the model in 
the commercial codes MSC.Nastran and ABAQUS. The commercial implementations of the ECTE 
model will be validated in comparison to a research code for the analysis cases of post-buckling, 
random response, and shape control of SMAHC beams. Extensive additional detail is available on 
the actual analysis input files and solution procedures in the companion document in reference 14 
and electronic versions of all input files are available at the URL http://stabserv.larc.nasa.gov . 

2. MODEL FORMULATION 

The salient features of the ECTE constitutive model for SMA and SMAHC materials are described 
in this section. A more detailed description can be found elsewhere. 5 There are various composite 
material types, structural fabrication techniques, and embedding strategies to integrate SMA 
actuators within a composite structure, e.g., between layers, within layers, etc. Likewise, there are 
many micromechanical models for evaluating the effective composite properties and resulting 
constitutive equations for such hybrid composites. The case considered here is that of a laminated 
composite composed of unidirectional pre-impregnated (pre-preg) layers with SMA actuators 
embedded within the laminae. A representative volume element from an individual SMAHC 
lamina, with principal material coordinates 1 and 2, is shown schematically in Figure 1 . This 
example is used for simplicity of illustration and it is understood that a similar approach could be 
used to model different composite materials, embedding strategies, and structure types. 
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The thermoelastic constitutive relations in principal material coordinates for a thin orthotropic 
lamina under conditions of plane stress can be shown to have the following form 
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where [Q\ is the reduced stiffness matrix, {s} is the total strain vector, and {a} is the vector of 
CTEs. The reduced stiffness terms are related to the effective engineering properties through the 
following equations. 
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Estimates of the effective engineering properties for the lamina in Figure 1 must be determined in 
order to make use of these constitutive relations. A mechanics of materials approach can be used as 
follows to develop the effective properties for such a SMAHC material. 


The most fundamental feature of the ECTE model is the axial constitutive relation for a SMA 
actuator, in which the non-mechanical (i.e., thennal and transformation) strain is represented by 
effective thermal strain: 


a (T) = E a{ T ) C -\l a \aW) dT 


(3) 


where E a (T ) is the Young’s modulus of the SMA material, s x is the total axial strain, a Ui is the 

effective coefficient of thennal expansion, and the subscript a indicates that the quantity is specific 
to the actuator material. It can be seen in Equation (3) that the ECTE model is intended to model 
thennally induced transformation phenomena, as opposed to stress induced (pseudoelasticity), which 
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is the restriction mentioned previously. It can be shown that the thennal strain in Equation (3) is 
governed by the usual thermoelastic effects at temperatures below the austenite start temperature 
( T < A s ) and is related to the recovery stress ( cr. ) and modulus (E a ) at higher temperatures 

( T>A s ): 

a r (T) = -E a (T)^a Wi (T)dT or a Ui (z)dr = -cr r (T) / E a (T) (4) 

This expression implies that measurement of the fundamental engineering properties ( a la for 
T < A s , cr. and E a for T > A s ) from experiments approximating the intended use condition (e.g., 

prestrain level and constraint) is all that is required to quantify the model. For example, an actuator 
under stiff constraint and restricted to small prestrain recovery in service would be well represented 
by measurements of blocked recovery stress. A similar constitutive expression can be written for the 
actuator in the transverse direction: 

(T 2a( T ) = E a( T ) £ 2~\ T r a 2a( T ) dT ( 5 ) 

where a 2a is not related to recovery stress, but is still nonlinear due to differing martensitic and 

austenitic properties. Note that it has been assumed that the Young’s modulus for the SMA is 
directionally independent. The analogous expressions for the composite matrix material, denoted by 
the subscript m, follow directly. 

°i,« = E lm CO - \ T ' a u„ (?)dr and a 2rn = E lm ( T ) a 2m (r)dr (6) 

It can be shown that one form for the effective properties of a SMAHC lamina can be constructed 
from these constituent relations as follows. 
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These effective engineering properties can be used in conjunction with Equations (1) and (2) to 
define the thermoelastic characteristics of a SMAHC lamina. 


Note that the ECTE model is consistent with classical thermoelasticity, so the model can be readily 
implemented in any finite element code that has the capability for analyzing laminated composite 
materials with temperature dependent material properties. It is envisioned that a similar approach 
might be possible for incorporation of a phenomenological or micromechanics-based 1-D 
constitutive model for the SMA in a SMAHC material model by replacing Equation (3) with the 
corresponding constitutive and kinetic relations. 

3. COMMERCIAL CODE IMPLEMENTATION 

One approach for deriving the constitutive relations governing the thermoelastic properties of a 
SMAHC lamina was shown in the previous section. Attention will now turn to methods for 
implementing material models of this type using the commercial pre-/post-processor MSC.Patran 
and generating the data structure necessary to perform various types of analyses on SMAHC 
structures within the framework of the commercial finite element codes MSC.Nastran and 
ABAQUS. 


3.1. MSC.Patran 
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Note in Equations (7) that the volume fraction of SMA for a given SMAHC layer can range from 
0% to 100%. This allows for a large variety of material systems and material system modeling 
approaches. For example, consider the case of a SMAHC laminate with SMA ribbon embedded 
within or between specific laminae. Furthermore, let the actuator cross-section aspect ratio be large; 
thickness on the order of the lamina thickness and width on the order of typical finite element (e.g., 
plate or shell element) dimensions. A typical finite element cross section with this material system 
type is shown schematically in Figure 2a, which indicates that individual finite elements may consist 
of 100% SMA and/or 100% host composite matrix layers. Note that unrealistically few layers are 
shown for clarity of illustration. This type of SMAHC element requires only the distinct properties 
of the constituent materials to define the properties of the individual layers in a particular finite 
element and will be referred to as the mon element (or laminate) type, indicative of the monolithic 
SMA inclusion within a layer. Evaluation of element properties with a material system of this type 
can be accomplished in MSC.Patran by defining 2-D orthotropic materials to represent the SMA and 
composite matrix material separately. Orthotropic engineering properties are typically known or 
measured for a particular composite matrix material. The full complement of orthotropic 
engineering properties for the SMA material can be constructed from measurements of Young’s 
modulus E a and recovery stress cr, by assuming that the Young’s modulus and Poisson’s ratio are 

directionally independent, the shear modulus is related to the Young’s modulus through the typical 
isotropic relation G = E / 2(1 + v ) , the axial CTE is detennined from the thennal strain in Equation 
(4), and the transverse CTE varies from martensitic to austenitic values. Finally, Classical 
Lamination Theory can be used to develop the SMAHC laminate stiffness coefficients, thermal 
forces, and thermal moments. 
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Alternatively, a similar structural functionality and perfonnance could be achieved with SMA 
actuators in the form of round wire with a diameter on the order of the lamina thickness, which 
suggests several actuators to a typical element width. A schematic of a finite element cross section 
with a material system of this type is shown in Figure 2b, where again relatively few layers are 
shown for clarity. It can be seen that individual finite elements may consist of all SMAHC layers 
with varying orientations and SMA content, combinations of traditional and SMAHC layers, etc. 
Thus, the material model for individual layers in a finite element of this type requires properties for 
mixtures of the constituents, defined by relations like those in Equations (7), and will be referred to 
as the mix element (or laminate) type, indicative of the mixture of SMA and host within a layer. In 
this case, evaluation of element properties can be accomplished in MSC.Patran by defining 3-D 
orthotropic materials representing the SMA and composite matrix material, each with equivalent 
mechanical properties in the two transverse directions. The constituent materials can then be 
combined to form a SMAHC material by the Rule-of-Mixtures approach. Lamination sequences of 
traditional composite and/or SMAHC laminae can be modeled through Classical Lamination Theory 
to develop the corresponding laminate stiffness coefficients, thermal forces, and thermal moments. 

While MSC.Patran supports output of temperature dependent material property data for isotropic 
and anisotropic materials, it does not for ortho tropic materials. Thus, modeling strategies for an 
anisotropic SMAHC material may be plausible for MSC.Patran output. The present study is 
primarily concerned with laminated SMAHC materials, which are composed of ortho tropic layers. 
MSC.Patran will produce the laminate definition information (PCOMP bulk data entries in 
MSC.Nastran and *SHELL SECTION options in ABAQUS) and the temperature independent 
material property data for the individual laminae (MATi bulk data entries in MSC.Nastran and 
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combinations of temperature independent versions of the 

* MATER I A L/* ELASTIC/* EX PANSION/ctc. options in ABAQUS) for both laminate types (man 
and mix ) described above. However, MSC.Patran will not produce the required temperature 
dependent material property data (MATTi and TABLEMi bulk data entries in MSC.Nastran and 
appropriate combinations of temperature dependent versions of the 

* MATER I A L/* ELASTIC/* EX PANSION/ctc. options in ABAQUS) for either laminate type. This 
limitation is not particularly consequential for elements of the mon laminate type, layers alternating 
between 100% SMA and 100% host composite, because a single tabulated data set for each 
constituent is all that is required. However, some additional effort is required for elements of the 
mix laminate type, where any given layer may consist of mixtures of SMA and host composite. In 
this case, a new tabulated dataset must be created from the constituent properties, by using 
Equations (7), for every material combination employed. The analysis input file must be hand 
altered, for either laminate type, to incorporate the appropriate temperature dependent material 
properties. Of course, this limitation might be removed in some other pre-processor that is unknown 
to the authors. 

3.2. MSC.Nastran 

The capability for analyzing laminated composite materials with temperature dependent material 
properties was recently added to MSC.Nastran. A brief description of this capability is given here, 
but further detail can be found in the MSC.Nastran Version 2005 Release Notes and Quick 
Reference Guide. This capability is available for QUAD4 and TRIA3 h-element types and is 
implemented for SOL 106, Nonlinear Static (NLSTATIC) solution procedure. Other response 
characteristics of SMAHC structures can be determined by a RESTART from any particular 
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increment in a nonlinear static solution sequence. This will be described in more detail in later 
sections. 

Two new parameters have been introduced to allow user specification of options involved in 
analysis of laminated composites with temperature dependent material properties. The parameter 
COMPMATT indicates the intent for analysis of laminated composites with temperature dependent 
properties and the parameter EPSILONT indicates the form of the thermal expansion property data. 
The form for specification of these parameters is as follows: PARAM, COMPMATT, YES (Default 
NO) and PARAM, EPSILONT, INTEGRAL (Default SECANT). Setting the parameter 
COMPMATT to YES allows PCOMP bulk data entries, which define the composite material 
laminate, to reference temperature dependent material properties in TABLEMi bulk data entries 
through corresponding MATTi (i= 1 ,2,8) and MATi (i=l,2,8) bulk data entries. The default value for 
the parameter EPSILONT causes MSC.Nastran to expect the secant form for CTE data in the 
TABLEMi bulk data entries and thennal strain is calculated according to the equation 

s(T) = a(T)(T - T ef )-a(T i )(T i - T ef ) (8) 

where the subscripts i and revindicate initial and reference quantities, respectively. In the event that 
the load temperature T is not equal to a tabulated value, a(T) is estimated by linear interpolation 
within the tabulated data. A value of INTEGRAL for the parameter EPSILONT indicates that the 
CTE data is in tangent form and the thennal strain should be calculated as the integral of the 
interpolated tangent CTE. Thus, the thermal strain developed over the temperature change 7j — » 71, 

within the tabulated data interval ( a j ,T i ) and (. 0 Cj,Tj ) is evaluated as 
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( 9 ) 



and the total thermal strain is the sum of integrals over each interval from the initial to the load 
temperature. Note that the thermal strain over the entire interval above yields 

t j { 

]a(T)dT = -(a J +a i )(T J -T i ) (10) 

1 } ^ 

In general, the integral approach has greater accuracy with increasing temperature interval size. 

The user also has the option of specifying thermal strain directly. This is indicated in the bulk data 
file by assigning a negative table identification number (TID) to the TABLEMi material property 
table for the corresponding thennal expansion data. Of course, the same negative TID must also 
appear in the MATTi bulk data entry referring to the temperature dependent thermal strain data. 

This capability may be convenient in many cases where thermal strain is known, which is 
particularly useful for the ECTE model because thennal strain in the actuator’s axial direction is 
known from Equation (4) for T > A s . In this case, temperature difference evaluation of thennal 

strain, i.e, by the secant or integral approach, would be necessary only if analysis were required at 
temperatures less than A s . 

3.3. ABAQUS 

ABAQUS is also capable of modeling laminated composite materials with temperature dependent 
material properties. This capability is available for many element types and solution procedures, so 
the user is referred to the ABAQUS User’s and Reference Manuals. A large variety of solution 
combinations can be achieved by making use of restart capabilities with the * RESTART option 
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and/or multiple analysis steps (*STEP option) within a single analysis input file. The *SHELL 
SECTION, COMPOSITE option is used to describe the material and orientation stacking sequence, 
which refers to corresponding *MATERIAL options. *ELASTIC, *EXPANSION, and *DENSITY 
options are then associated with each * MATERIAL option to fully prescribe the temperature 
dependent material properties. ABAQUS accepts thermal expansion data in terms of secant CTE 
values, but not in terms of tangent CTE values. ABAQUS also accepts input of thennal strain data 
by employing the USER parameter on the *EXPANSION option and defining the user subroutine 
UEXPAN, but this capability is not explored in this study. 


4. EXAMPLE PROBLEMS 

The commercial implementations of the ECTE model will now be demonstrated and validated 
through comparison with a research code for two relatively simple SMAHC structures. The cases to 
be considered include a beam clamped at both ends and a cantilever beam. Both mon and mix 
laminate types will be considered in the subsequent analyses of static and dynamic structural 
response. 

4.1. Clamped SMAHC Beam 

Consider a beam with a length and width of 45.72x2.54 cm (18x1 inches) and clamped boundary 
conditions at both ends. The lamination stacking sequence is (45/0/-45/90)2 S with SMA ribbon 
material replacing a width of 1.143 cm (0.45 inches) of the 0° glass-epoxy layers about the beam 
centerline, bounded by 0.6985 cm (0.275 inch) wide strips of glass-epoxy on either side. The 
thickness of each glass-epoxy layer is 0.012383 cm (0.004875 inches) and that for the SMA is 
0.01524 cm (0.006 inches). The motivation for this particular geometry and material system is from 
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a previous experimental validation study of the ECTE model as implemented in a research code. 

The finite element mesh used to model the beam specimen consists of 36x4 rectangular elements 
with two element widths modeling the SMA inclusion and one element width modeling each of the 
glass-epoxy strips along the outer edges of the beam. A graphical representation of the beam 
specimen and finite element mesh is shown in Figure 3. QUAD4 quadrilateral shell elements were 
used in MSC.Nastran, S4 quadrilateral shell elements were used in ABAQUS, and rectangular plate 
elements 15 were used in the research code. Note that a mon laminate type modeling approach is used 
in this case, which results in two element types in the finite element model; one SMAHC and one 
glass-epoxy. The two element types have differing overall thicknesses, 0.198 cm (0.078 inches) for 
the glass-epoxy and 0.210 cm (0.0825 inches) for the SMAHC, but the same mid-plane. The 
temperature dependent material properties for the two constituents (Nitinol and glass-epoxy) are 
given in Table 1 and Table 2, where the thermal expansion data is given in secant CTE form. 

The physical conditions that are to be numerically simulated for this specimen involve imposing a 
uniformly distributed thermal load of 121.1°C (250°F) in the presence of the clamped boundary 
conditions and subjecting the beam to random out-of-plane inertial load. The inertial load is 
characterized by a 0.25g RMS acceleration uniformly distributed over a bandwidth of 0-400 Hz, i.e., 
band-limited white noise with a spectral level of 150.51 (cm/ s") /Hz (23.33 (in/s") /Hz). It is clear 
that the combination of thennal loads and immovable in-plane boundaries presents the opportunity 
for thermoelastic instability and large thermal post-buckling deflections. Typical numerical 
approaches for nonlinear static solutions involving stability require some seeding of the solution by 
an initial deflection to prevent numerical ill conditioning. This initial deflection can be constructed 
in many ways including weighted combinations of buckling eigenvectors, geometric imperfections, 
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and other static deflections due to real or fictitious loads. The approach taken here is that of 
introducing a static gravity load of lg (9.81 m/s 2 , 386.4 in/s 2 ), nonnal to the plane of the beam, in a 
separate static analysis and using the resulting deflections as geometric imperfections. 

The required solutions were achieved within the framework of MSC.Nastran by perfonning a 
nonlinear static solution (SOL 106) with a uniform thermal load increasing from ambient 
temperature to 121.1°C (250°F). The geometric imperfections (gravity deflections) were introduced 
into the nonlinear static solution by perturbing the coordinates (i.e., x+u, y+v, z+w) of the GRID 
bulk data entries defining the nodal locations. Linear random analyses were performed, subsequent 
to the nonlinear static solution, using a modal approach incorporating the first 10 modes, a critical 
damping ratio of 0.5% for all modes, and a frequency resolution of 0.25 Hz at four temperatures; 
23.9°C (75 °F), 65.6°C (150°F), 93.3°C (200°F), and 121. 1°C (250°F). Linear random analysis was 
accomplished as a variant of frequency response analysis (SOL 111) with the definition of the 
RANDPS bulk data entry that referred to a frequency dependent spectral load given by a TABRND1 
bulk data entry. The imperfections were again introduced into the MSC.Nastran reference 
temperature, random input file by perturbing the nodal coordinates. Linear random analysis at 
elevated temperature was accomplished by using a RESTART statement to reuse the database from 
the previous nonlinear static analysis and by specifying the specific thermal load increment number 
with the parameter NMLOOP, i.e., with the bulk data entry PARAM,NMLOOP,#. Thus, the 
geometric imperfections and thennoelastic equilibrium state are implicitly included in the elevated 
temperature random analyses because of the restart from the nonlinear static solution. 
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A similar analysis sequence was achieved by a slightly different approach in ABAQUS. A single 
analysis input file was created and several analysis steps were defined to perfonn the analysis 
sequentially. The nonlinear static solution was thereby progressed through the desired equilibrium 
conditions and the random response was determined as linear perturbation steps at the specified 
thennoelastic equilibrium states, i.e. 23.9°C (75°F), 65.6°C (150°F), 93.3°C (200°F), 121. 1°C 
(250°F). The geometric imperfections were introduced in the ABAQUS input file by using the 
* IMPERFECTION, FILE option to read the deflection data from the gravity analysis results file. 

All of the linear random response steps were again performed using the modal approach 
incorporating the first 10 modes, each with 0.5% critical damping ratio, and with a frequency 
resolution of 0.25 Hz. 

A comparison of the nonlinear static analysis results from MSC.Nastran, ABAQUS, and the research 
code is shown in Figure 4, which shows the mid-span out-of-plane deflection as a function of 
temperature using secant CTE thermal strain calculation. It can be seen that excellent agreement has 
been achieved between all three implementations of the ECTE model. Analogous results are shown 
in Figure 5 for two alternate thennal strain calculation cases. In the integral case, thermal strain is 
taken directly from tabulated data for the axial direction of the SMA (Equation (4)) and calculated 
by the integral method with tangent CTE data otherwise. Results for the strain case are generated by 
using thennal strain data for both directions in both materials; SMA and glass-epoxy. It can be seen 
that the results are essentially equivalent to the secant CTE results and excellent agreement between 
the codes is again achieved. Note that analysis was not possible in ABAQUS for the integral case 
and possible but not performed for the strain case. The results in both figures show that traditional 
thennoelastic behavior prevails up to a temperature of approximately 37.8°C (100°F), where the 
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SMA actuation authority begins to dominate the stress state. At higher temperatures, the structure is 
rendered flat, even eliminating the initial imperfection (gravity) deflections. The stiffening effect 
dominates over a large temperature range, implying a significant effect on the dynamic response. 
More detailed descriptions of this complex thermoelastic behavior can be found in reference 5. 

Predictions from the three codes for the beam’s linear random response at the two temperature 
extremes are shown in Figure 6 and Figure 7. The figures show the mid-span displacement power 
spectral density (PSD) versus frequency. Again, excellent correlation between the three codes is 
observed, as expected from the agreement between the nonlinear static analysis results. These 
results clearly demonstrate the large stiffening effect afforded by the SMA, which results in large 
shifts in the resonant frequencies and reductions in the dynamic response amplitude. For example, 
the fundamental frequency and RMS displacement response at 121.1°C (250°F) are 4.37/9.49 times 
higher/lower than their reference temperature counterparts. Additional insight into the dynamic 
behavior of SMAHC structures is given in reference 13. 

Similar analyses were conducted on a SMAHC beam specimen of a mix laminate type. Consider the 
same specimen geometry and layer orientations, but with a unifonnly distributed mixture of SMA 
and glass-epoxy in the 0° layers. This configuration is intended to simulate inclusion of small 
diameter round wire in the composite laminate. In this case, all layers were taken to have a 
thickness of 0.012383 cm (0.004875 inches), which gives an overall thickness of 0.198 cm (0.078 
inches) for all elements. A volume fraction of 55.38% SMA was specified for the 0° layers to 
render an overall SMA volume fraction equivalent to the previous mon laminate case. Equations (7) 
were used with the specified volume fractions and the properties in Table 1 and Table 2 to generate 
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effective properties for the 0° SMAHC layers, shown in Table 3. These SMAHC effective 
properties were used, along with those for glass-epoxy alone, to construct the laminate properties for 
the model in Figure 3. Note that all elements in this finite element model share the same SMAHC 
properties of the mix laminate type, in contrast to the previous beam with SMA ribbon. A 
comparison of the thermal post-buckling (nonlinear static) results for this case is shown in Figure 8, 
where it can be seen that the maximum deflection has increased marginally and similar excellent 
agreement between the codes has been achieved. Similar performance between the two SMAHC 
laminate types was achieved intentionally in this case, but the utility of the modeling approach in 
achieving optimal designs is implied. 

4.2. Cantilever SMAHC Beam 

Focus is now directed to shape/deflection control of a cantilever SMAHC beam. Consider a 
cantilever SMAHC beam with a length and width of 22.86x2.54 cm (9x 1 inches), or half of the 
model in Figure 3. A mon laminate case similar to that discussed previously is taken with the same 
stacking sequence (45/0/-45/90)2 S , but with SMA ribbon material replacing a width of 1.143 cm 
(0.45 inch) about the beam centerline in only the 2 nd 0° glass-epoxy layer. Again, the outer 0.699 
cm (0.275 inch) strips of the beam are modeled as solely glass-epoxy and the thickness of the glass- 
epoxy and SMA are taken to be 0.012383 cm (0.004875 inch) and 0.015 cm (0.006 inch), 
respectively. Two element types are again required, one SMAHC and one glass-epoxy only, with 
overall thicknesses of 0.201 cm (0.079 inches) and 0.198 cm (0.078 inches), respectively. The 
material properties for the laminate constituents are given in Table 1 and Table 2. The finite element 
model is subjected to a uniformly distributed thermal load of 121.1°C (250°F) to study the deflection 
response. No seeding of the deflection is necessary in this case due to the development of a thermal 
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moment. Excellent agreement between MSC.Nastran and ABAQUS is apparent in the plot of tip 
displacement versus temperature shown in Figure 10. Excellent agreement between the three 
thennal strain calculation methods in MSC.Nastran is also demonstrated. 

A similar comparison is shown in Figure 1 1 for the analogous mix laminate case of uniformly 
distributed SMA throughout the 2 nd 0° layer, which was modeled by a 55.38% SMA volume fraction 
in the single SMAHC layer of 0.012383 cm (0.004875 inch) thickness. In this case all finite 
elements have the same lamination, overall thickness of 0.198 cm (0.078 inches), and material 
properties given in Table 1 and Table 3. It can be seen that excellent correlation is achieved 
between MSC.Nastran and ABAQUS for this case also. These results again demonstrate the 
actuation authority of the SMA and indicate that large deflections are possible using relatively little 
SMA. Results from the research code are excluded from the deflection control comparisons because 
the response is beyond the scope of the formulation. 

The deflection control example provides an introduction to the concept of shape control of structures 
by embedded SMA actuation. This topic is currently of significant interest for a variety of advanced 
aerospace vehicle applications, e.g., see references 16 and 17. Note that use of the ECTE model for 
shape control applications should be limited to cases in which the recovered strain of the SMA 
actuators is small relative to the prestrain, unless accommodations are made to schedule diminishing 
actuation authority. Also, prediction of the response during unloading requires data in addition to 
that in Table 1-Table 3 because of the hysteresis inherent to SMA materials. 
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Extensive additional detail on the analysis input files and solution procedures for the cases presented 
in this study are available in reference 14. ASCII listings of all of the input files and data necessary 
to duplicate the results in this document are also available in that reference, as well as in electronoic 
form at the URL http://stabserv.larc.nasa.gov . 


5. SUMMARY AND CONCLUSIONS 

A thermoelastic model for shape memory alloy (SMA) and SMA hybrid composite (SMAHC) 
materials has been implemented in the commercial finite element codes MSC.Nastran and 
ABAQUS. The model is based upon definition of an effective coefficient of thermal expansion 
(ECTE). The model can be readily implemented in any code that has capability for structural 
analysis with temperature dependent material properties. The model was applied to two different 
types of SMAHC laminated materials, corresponding to inclusion of SMA actuators in ribbon and 
round wire product forms. The SMA ribbon inclusion material model resulted in finite elements 
with layers composed of either host composite material or SMA material alone. The model for 
inclusion of round SMA wire required calculation of effective properties for mixtures of host 
composite and SMA in the finite element layers. This calculation can be performed automatically in 
a suitable finite element pre-processor. Alternatively, the effective properties can be calculated 
externally, e.g., using the formulation presented in this document, from the constituent properties 
and hand entered into the analysis input files. 

SMAHC structural analysis was demonstrated on SMAHC beams of two types; a beam clamped at 
both ends and a cantilever beam. Results from nonlinear static (post-buckling) and random vibration 
responses of the first specimen type showed excellent correlation between the two commercial 
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implementations and the research code. Complex thermoelastic behavior and the enormous control 
authority of the SMA actuators were demonstrated. Three different forms of thermal strain 
calculation in MSC.Nastran were shown to produce equivalent results, as expected. Analysis of the 
two different SMAHC material configurations, one involving large cross section aspect ratio SMA 
ribbon and another involving small diameter SMA wire, showed similar response and equally good 
analysis comparison with implications for parametric studies and/or optimal design. Analysis of the 
cantilever specimen introduced the use of the ECTE model for shape control prediction. Again, 
results are presented from analysis of two material configurations. The deflection of the cantilever 
beam in response to heating demonstrated the efficient control authority of the embedded SMA. 
Excellent correlation was achieved between the two commercial codes. 
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SYMBOLS 


A Transformation temperature 

E Elastic modulus 

G hear modulus 

Q Composite lamina reduced stiffness 

T Temperature 

Greek 

a Coefficient of thermal expansion 

y Engineering shear strain 

s Normal strain 

v Poisson’s ratio 

a Normal stress 

x Shear stress 

Subscripts 

1,2,6 Orthotropic tensor indices in contracted notation 

a Associated with actuator material 

m Associated with matrix material 

o Ambient or reference quantity 

ref Reference quantity 

s Indicates transformation start 
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Figure 1: Representative volume element for a SMAHC lamina. 

Figure 2: Representative beam cross sections for the mon (a) and mix (b) SMAHC laminate types. 
Figure 3: SMAHC beam clamped at both ends with finite element mesh. 

Figure 4: Post-buckling deflection versus temperature for the mon laminate case using secant CTE 
data. 

Figure 5: Post-buckling deflection versus temperature for the mon laminate case using combination 
tangent CTE/thermal strain data and all-thermal strain data. 

Figure 6: Random inertial mid-span displacement response PSD at 23.9°C (75°F) for the mon 
laminate case. 

Figure 7: Random inertial mid-span displacement response PSD at 121.1°C (250°F) for the mon 
laminate case. 

Figure 8: Post-buckling deflection versus temperature for the mix laminate case using secant CTE 
data. 

Figure 9: SMAHC cantilever beam with finite element mesh. 

Figure 10: Tip deflection versus temperature for the mon laminate case using secant CTE, 
combination tangent CTE/thermal strain data, and all-thermal strain data types. 

Figure 1 1 : Tip deflection versus temperature for the mix laminate case using secant CTE, 
combination tangent CTE/thennal strain, and all-thermal strain data types. 
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Figure 1: Representative volume element for a SMAHC lamina 
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Figure 2: Representative beam cross sections for the mon (a) and mix (b) SMAHC laminate types. 


Host 
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Figure 3: SMAHC beam clamped at both ends with finite element mesh. 
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Mid-span Deflection, cm 



Figure 4: Post-buckling deflection versus temperature for the man laminate case using secant CTE data. 
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Mid-span Deflection, cm 



Figure 5: Post-buckling deflection versus temperature for the mon laminate case using combination tangent CTE/thermal 

strain data and all-thermal strain data. 
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Mid-span Deflection, cm 



Figure 8: Post-buckling deflection versus temperature for the mix laminate case using secant CTE data. 
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Figure 9: SMAHC cantilever beam with finite element mesh. 
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Tip Deflection, cm 



Figure 10: Tip deflection versus temperature for the mon laminate case using secant CTE, combination tangent 

CTE/thermal strain, and all-thermal strain data types. 
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Tip Deflection, cm 



Figure 11: Tip deflection versus temperature for the mix laminate case using secant CTE, combination tangent 

CTE/thermal strain, and all-thermal strain data types. 
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Table 1: Temperature dependent orthotropic material properties for glass-epoxy (p=5.259E-3 g/cm 3 , secant CTE data 

presented with T re f=23.9°C). 


T, °C 

a-,, l°C 

a 2 , 1° C 

ro 

Q- 

uT 

E 2 , Pa 

Gi2< Pa 

V 12 

15.6 

5.270E-06 

1.105E-05 

4.930E+10 

2.000E+10 

9.653E+09 

0.29 

21.1 

5.373E-06 

1.155E-05 

4.930E+10 

2.000E+10 

9.653E+09 

0.29 

26.7 

5.679E-06 

1.306E-05 

4.930E+10 

2.000E+10 

9.653E+09 

0.29 

32.2 

5.963E-06 

1.480E-05 

4.923E+10 

1.972E+10 

9.446E+09 

0.29 

37.8 

6.248E-06 

1.654E-05 

4.916E+10 

1.944E+10 

9.239E+09 

0.29 

43.3 

6.433E-06 

1.788E-05 

4.909E+10 

1.920E+10 

9.067E+09 

0.29 

48.9 

6.619E-06 

1 .922E-05 

4.902E+10 

1.896E+10 

8.895E+09 

0.29 

54.4 

6.694E-06 

2.003E-05 

4.892E+10 

1.872E+10 

8.722E+09 

0.29 

60.0 

6.770E-06 

2.083E-05 

4.882E+10 

1.848E+10 

8.550E+09 

0.29 

65.6 

6.788E-06 

2.131E-05 

4.875E+10 

1.820E+10 

8.412E+09 

0.29 

71.1 

6.779E-06 

2.164E-05 

4.875E+10 

1.779E+10 

8.274E+09 

0.29 

76.7 

6.751 E-06 

2.180E-05 

4.871 E+10 

1 .741 E+10 

8.102E+09 

0.29 

82.2 

6.723E-06 

2.196E-05 

4.868E+10 

1.703E+10 

7.929E+09 

0.29 

87.8 

6.688E-06 

2.200E-05 

4.864E+10 

1.662E+10 

7.757E+09 

0.29 

93.3 

6.653E-06 

2.203E-05 

4.861 E+10 

1.620E+10 

7.585E+09 

0.29 

98.9 

6.630E-06 

2.202E-05 

4.861 E+10 

1.576E+10 

7.171 E+09 

0.29 

104.4 

6.608E-06 

2.201 E-05 

4.861 E+10 

1 .531 E+10 

6.757E+09 

0.29 

110.0 

6.606E-06 

2.203E-05 

4.858E+10 

1.486E+10 

6.378E+09 

0.29 

115.6 

6.604E-06 

2.205E-05 

4.854E+10 

1 .441 E+10 

5.999E+09 

0.29 

121.1 

6.619E-06 

2.210E-05 

4.854E+10 

1.400E+10 

5.585E+09 

0.29 

126.7 

6.644E-06 

2.219E-05 

4.861 E+10 

1.345E+10 

5.171 E+09 

0.29 

132.2 

6.675E-06 

2.232E-05 

4.864E+10 

1.293E+10 

4.723E+09 

0.29 

137.8 

6.707E-06 

2.245E-05 

4.868E+10 

1 .241 E+10 

4.275E+09 

0.29 

143.3 

6.744E-06 

2.260E-05 

4.875E+10 

1.189E+10 

3.861 E+09 

0.29 

148.9 

6.781 E-06 

2.275E-05 

4.882E+10 

1.138E+10 

3.448E+09 

0.29 
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Table 2: Temperature dependent orthotropic material properties for Nitinol (p=1.481E-2 g/cnT, secant CTE data 

presented with T re f=23.9°C). 

T, °C a^l°C a 2 ,l°C E 1( Pa E 2 , Pa G 12 , Pa v 12 

15.6 6.606E-06 6.606E-06 2.717E+10 2.717E+10 1.045E+10 0.3 

21.1 6.606E-06 6.606E-06 2.717E+10 2.717E+10 1.045E+10 0.3 

26.7 -3.089E-05 6.606E-06 2.480E+10 2.480E+10 9.538E+09 0.3 

32.2 -3.871 E-05 6.606E-06 2.243E+10 2.243E+10 8.628E+09 0.3 

37.8 -5.040E-05 6.606E-06 2.006E+10 2.006E+10 7.717E+09 0.3 

43.3 -5.705E-05 6.606E-06 2.570E+10 2.570E+10 9.886E+09 0.3 

48.9 -1.005E-04 6.676E-06 3.134E+10 3.134E+10 1.206E+10 0.3 

54.4 -1.345E-04 6.834E-06 3.698E+10 3.698E+10 1.423E+10 0.3 

60.0 -1.416E-04 7.040E-06 4.262E+10 4.262E+10 1.639E+10 0.3 

65.6 -1.340E-04 7.275E-06 4.827E+10 4.827E+10 1.856E+10 0.3 

71.1 -1.196E-04 7.529E-06 5.486E+10 5.486E+10 2.110E+10 0.3 

76.7 -1.053E-04 7.795E-06 6.145E+10 6.145E+10 2.363E+10 0.3 

82.2 -9.932E-05 8.070E-06 6.421E+10 6.421E+10 2.469E+10 0.3 

87.8 -9.898E-05 8.325E-06 6.313E+10 6.313E+10 2.428E+10 0.3 

93.3 -9.719E-05 8.538E-06 6.206E+10 6.206E+10 2.387E+10 0.3 

98.9 -9.127E-05 8.721E-06 6.392E+10 6.392E+10 2.458E+10 0.3 

104.4 -8.458E-05 8.878E-06 6.578E+10 6.578E+10 2.530E+10 0.3 

110.0 -7.877E-05 9.015E-06 6.764E+10 6.764E+10 2.602E+10 0.3 

115.6 -7.335E-05 9.135E-06 6.950E+10 6.950E+10 2.673E+10 0.3 

121.1 -6.816E-05 9.241 E-06 7.136E+10 7.136E+10 2.745E+10 0.3 

126.7 -6.576E-05 9.336E-06 7.084E+10 7.084E+10 2.725E+10 0.3 

132.2 -6.325E-05 9.421E-06 7.032E+10 7.032E+10 2.704E+10 0.3 

137.8 -6.104E-05 9.498E-06 6.979E+10 6.979E+10 2.684E+10 0.3 

143.3 -5.890E-05 9.568E-06 6.927E+10 6.927E+10 2.664E+10 0.3 

148.9 -5.71 1 E-05 9.632E-06 6.874E+10 6.874E+10 2.644E+10 0.3 
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Table 3: Temperature dependent orthotropic material properties for SMAHC mixture with 55.38% SMA content 
(p=1.055E-2 g/cm 3 , secant CTE data presented with T re f=23.9°C). 

T, °C a-i, /°C a 2 , /°C E 1( Pa E 2 , Pa G 12 , Pa v 12 

15.6 5.813E-06 8.589E-06 3.704E+10 2.342E+10 1.008E+10 0.296 

21.1 5.874E-06 8.812E-06 3.704E+10 2.342E+10 1.008E+10 0.296 

26.7 -8.380E-06 9.483E-06 3.573E+10 2.240E+10 9.589E+09 0.296 

32.2 -1.018E-05 1.026E-05 3.439E+10 2.114E+10 8.975E+09 0.296 

37.8 -1.280E-05 1.104E-05 3.305E+10 1.978E+10 8.329E+09 0.296 

43.3 -1.858E-05 1.164E-05 3.614E+10 2.233E+10 9.503E+09 0.296 

48.9 -4.077E-05 1.227E-05 3.923E+10 2.427E+10 1.041E+10 0.296 

54.4 -6.164E-05 1.272E-05 4.231 E+10 2.577E+10 1.110E+10 0.296 

60.0 -7.039E-05 1.319E-05 4.539E+10 2.693E+10 1.163E+10 0.296 

65.6 -7.082E-05 1.354E-05 4.848E+10 2.779E+10 1.207E+10 0.296 

71.1 -6.688E-05 1.382E-05 5.213E+10 2.843E+10 1.247E+10 0.296 

76.7 -6.160E-05 1.404E-05 5.577E+10 2.887E+10 1.274E+10 0.296 

82.2 -5.91 IE-05 1.427E-05 5.728E+10 2.872E+10 1.271 E+10 0.296 

87.8 -5.851 E-05 1.442E-05 5.667E+10 2.807E+10 1.245E+10 0.296 

93.3 -5.701 E-05 1.456E-05 5.606E+10 2.743E+10 1.219E+10 0.296 

98.9 -5.408E-05 1.466E-05 5.709E+10 2.704E+10 1.180E+10 0.296 

104.4 -5.055E-05 1.474E-05 5.812E+10 2.662E+10 1.137E+10 0.296 

110.0 -4.748E-05 1.482E-05 5.913E+10 2.617E+10 1.096E+10 0.296 

115.6 -4.456E-05 1.490E-05 6.015E+10 2.569E+10 1.052E+10 0.296 

121.1 -4.169E-05 1.498E-05 6.118E+10 2.523E+10 9.994E+09 0.296 

126.7 -3.999E-05 1.507E-05 6.092E+10 2.439E+10 9.381E+09 0.296 

132.2 -3.823E-05 1.518E-05 6.065E+10 2.359E+10 8.700E+09 0.296 

137.8 -3.667E-05 1.527E-05 6.037E+10 2.279E+10 8.000E+09 0.296 

143.3 -3.515E-05 1.538E-05 6.01 1 E+10 2.197E+10 7.335E+09 0.296 

148.9 -3.386E-05 1.549E-05 5.985E+10 2.115E+10 6.651E+09 0.296 
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